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Abstract 

This paper is concerned with the Wigner-Poisson-Fokker-Planck system, a kinetic evolu- 
tion equation for an open quantum system with a non-linear Hartree potential. Existence, 
uniqueness and regularity of global solutions to the Cauchy problem in 3 dimensions are es- 
tablished. The analysis is carried out in a weighted L 2 -space, such that the linear quantum 
Fokker-Planck operator generates a dissipative semigroup. The non-linear potential can be 
controled by using the parabolic regularization of the system. 

The main technical difficulty for establishing global-in-time solutions is to derive a-priori 
estimates on the electric field: Inspired by a strategy for the classical Vlasov-Fokker-Planck 
equation, we exploit dispersive effects of the free transport operator. As a "by-product" we 
also derive a new a-priori estimate on the field in the Wigner-Poisson equation. 
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1 Introduction 



The goal of this paper is to prove the existence and uniqueness of global-in-time solutions 
to the coupled Wigner-Poisson-Fokker-Planck (WPFP) system in three dimensions. This 
kinetic equation is the quantum mechanical analogue of the classical Vlasov-Poisson-Fokker- 
Planck (VPFP) system, which models the diffusive transport of charged particles (in plasmas 
e.g.). 

Wigner functions provide a kinetic description of quantum mechanics (cf. |Wij ) and have 
recently become a valuable modeling and simulation tool in fields like semiconductor device 
modeling (cf. |MRS| and references therein), quantum Brownian motion, and quantum optics 
( |CH iDi] ^. The real-valued Wigner function w(x,v,t) is a probabilistic quasi-distribution 
function in the position-velocity (x, v) phase space for the considered quantum system at 
time t. 

Its temporal evolution is governed by the Wigner-Fokker-Planck (WFP) equation: 

wt + v • V x w — @[V]w = f3div v (vw) + aA v w + 2jdW v (V x w) + aA x w, t > 0, (1-1) 
on the phase space x S IR 3 , v S IR 3 , with the initial condition 



With a vanishing right hand side, equation 1)1.1(1 would be the (diffusion-free) Wigner equa- 
tion. It describes the reversible evolution of a quantum system under the action of a (possibly 
time-dependent) electrostatic potential V = V(x,t). The potential effect enters in the equa- 
tion via the pseudo-differential operator Q[V]: 



where SV(x, rj, t) = V(x + U, t) — V(x — U,t) and T v ^w denotes the Fourier transform of w 
with respect to v: 



For simplicity of the notation we have set here the Planck constant, particle mass and charge 
equal to unity. 

The right hand side of <|1.1|) is a Fokker-Planck type model for the non-reversible interaction 
of this quantum system with an environment, e.g. the interaction of electrons with a phonon 
bath (cf. |CLI ICEFMj for derivations from reversible quantum systems, and |GGKS1 IStoj 
for applications in quantum transport). In 1)1.1)1 . (3 > is the friction parameter and the 
parameters a, 7 > 0, a > constitute the phase-space diffusion matrix of the system. In 
the kinetic Fokker-Planck equation of classical mechanics (cf. |Ril ICSVj ) one would have 
a = 7 = 0. For the WFP equation 1)1.1(1 we have to assume 



w(x, v, t = 0) 



w (x,v). 



(&[V]w)(x,v,t) 




i[V(t, x + — V„) - V(t, x - —V v )]w(x, v, t) 
1% li 

, n !o /2 / SV(x,r],t)J : ' v ^ v w(x,ri,t)e m ' v drj 
— ! I dV(x,ri,t)w(x,v,t)e i( - v - v '> v dvdri, 



(1.2) 
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which guarantees that the system is quantum mechanically correct. More precisely, it guar- 
antees that the corresponding von Neumann equation is in Lindblad form and that the 
density matrix of the quantum system stays a positive operator under temporal evolution 
(see jALMSl for details). In the sequel we shall therefore assume 

aa > 7 2 + — and aa > 7 2 . (1-3) 
16 

Hence, the principle part of the Fokker-Planck term is uniformly elliptic. This makes the 
present work complementary to |ALMSj . where the friction- free, hypoelliptic case (with 
a = f3 = 7 = 0) was analyzed. 

The WFP equation (jl.lj) is self-consistently coupled with the Poisson equation for the (real- 
valued) potential V = V[w](x,t): 

- AV = n[w], x£Wi 3 , t>0, (1.4) 

with the particle density 

„M(=M) := / «».,,«)*. (1.5) 
Jn 3 

This potential models the repulsive Coulomb interaction within the considered particle sys- 
tem in a mean-field description. 

The main analytical challenge for tackling Wigner-Poisson systems is the proper definition 
of n[w] in appropriate LP spaces. Due to the definition of the operator G in Fourier space, 
w G L 2 (JR 3 x IR 3 ) is the natural set-up. Without further assumptions, of course, this does 
not justify to define n[w]. We shall now summarize the existing literature of this field and 
the typical strategies to overcome the above problem: 

a) The standard approach for the Wigner-Poisson equation is to reformulate it as a Schro- 
dinger- Poisson system, where the particle density then appears in L 1 (cf. t BMj ICalj for the 
3D-whole space case). 

b) In one spatial dimension with periodic boundary conditions in x the Wigner-Poisson system 
(and WPFP) can be dealt with directly on the kinetic level. For w in a weighted L 2 -space, 
the nonlinear term 0[V]u; is then bounded and locally Lipschitz |AR| lACDj . The same 
strategy was also used in |Maj for the Wigner-Poisson system on a bounded (spatial) domain 
in three dimensions (local- in-time solution). 

c) By adapting L 1 -techniques from the classical Vlasov-Fokker-Planck equation, the 3D 
Wigner-Poisson-Fokker-Planck system was analyzed in |ALMSj (local-in-time solution for 
the friction-free problem) and |CLN| (global-in-time solution). The latter paper, however, is 
not a purely kinetic analysis as it requires to assume the positivity of the underlying density 
matrix. In both cases the dissipative structure of the system allows to control n[w]. 

d) In | Aril IAS| the 3D Wigner-Poisson and WPFP systems were reformulated as von- 
Neumann equations for the quantum mechanical density matrix. This implies n G L 1 (IR 3 ). 
While this approach is the most natural, both physically and in its mathematical structure, it 
is restricted to whole space cases. Extensions to initial-boundary value problems (as needed 
for practical applications and numerical analysis) seem unfeasible. 

e) For the classical Vlasov-Poisson-Fokker-Planck equation there exists a vast body of math- 
ematical literatur from the 1990's (cf. |Boll IBo2l ICS VI ICarl ICa2| ). and many of those tools 
will be closely related to the present work. 

In spite of the various existing well-posedness results for the WPFP problem, there is a need 
for a purely kinetic analysis, and this is our goal here. Such an approach could possibly allow 
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for an extension to boundary-value problems in the Wigner framework (where the positivity 
of the related density matrix is a touchy question). 

Mathematically we shall develop the following new tools and estimates that could be im- 
portant also for other quantum kinetic applications: In all of the existing literature on 
Wigner- Poisson problems (except jStej ) the potential V is bounded, which makes it easy to 
estimate the operator 0[V] in I?. Our framework for the local in time analysis does not 
yield a bounded potential. However, the operator only involves 5V, a potential difference, 
which has better decay properties at infinity. This observation gives rise to new estimates 
that are crucial for our local-in-time analysis. 

In order to establish global-in-time solutions we shall extend dispersive tools of Lions, 
Perthame and Castella (cf. |LP| IPe| ICPj for applications to classical kinetic equation) to 
the WP and WPFP systems. The fact that the Wigner function w also takes negative values 
gives rise to an important difference between classical and quantum kinetic problems: In the 
latter case, the conservation of mass and energy or pseudo-conformal laws do not provide 
useful a-priori estimates on w. We shall hence assume that the initial state lies in a weighted 
L 2 -space, but we shall not require that our system has finite mass or finite kinetic energy. 
Since the energy balance will not be used, this also implies that the sign of the interaction 
potential does not play a role in our analysis. 

This paper is organized as follows: In Section 2 we introduce a weighted L 2 -space for the 
Wigner function w that allows to define n[w] and the nonlinear term 0[V]io. In §3 we obtain 
a local-in-time, mild solution for WPFP using a fixed point argument and the parabolic 
regularization of the Fokker-Planck term. In §4 we establish a-priori estimates to obtain 
global-in-time solutions. The key point is to derive first L p -bounds for the electric field 
W by exploiting dispersive effects of the free kinetic transport. "Bootstraping" then yields 
estimates on the Wigner function in a weighted L 2 -space. Finally, we give regularity results 
on the solution. The technical proofs of several lemmata are defered to the Appendix. 
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2 The functional setting 

In this section we shall discuss the functional analytic preliminaries for studying the non- 
linear problem (|1.1|) - (|1.5|) . First we shall introduce an appropriate "state space" for the 
Wigner function w which allows to "control" the particle density n[w] and the self consistent 
potential V[tw]. Next, we shall discuss the linear Wigner- Fokker-Planck equation and the 
dissipativity of its (evolution) generator A. 

2.1 State space and selfconsistent potential 

Let us introduce the following weighted (real valued) L 2 -space 

X := L 2 (1R 6 ;(1 + \v\ 2 ) 2 dx dv) , (2.1) 
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endowed with the scalar product 

<u,w> x = / / uw{\ + \v\ 2 ) 2 dxdv. (2.2) 

The following proposition motivates the choice of X as the state space for our analysis. 

Proposition 2.1 For all w G X, the function n[w] defined by n[w](x) := fw(x,v) dv, 
x G IR 3 , belongs to L 2 (IR 3 ) and satisfies 

II^MIIl^r 3 ) < Clkllx, (2.3) 
with a constant C independent of w. 

Here and in the sequel C shall denote generic, but not necessarily equal, constants. 
Proof. By using Holder inequality in the v-integral, we get 



InHH^S) < /( /|tp(s,i;)| a (l + M a ) a (fo) ( / (1+ ^ |2)2 ) dx = C\\w\\ 2 x . 



□ 



Remark 2.2 The choice of the \v\ 2 weight was already seen to be convenient to control the 
L 2 -norm of the density on a bounded domain o/IR 3 (cf. IMcfl ). 

The subsequent analysis would hold also by including a symmetric weight in the x-variable 
(i.e. for w G L 2 (IR 6 ; (1 + \x\ 2 + \v\ 2 ) 2 dx dv) ), which would yield a L p -bound with p G (3/2,2] 
for the density. On the other hand, Lemmata \2.f\ \2.1CA would prevent us from introducing a 
non- symmetric weight in x. 

In this framework the following estimates for the self-consistent potential hold. 

Proposition 2.3 For all w S X, the (Newton potential) solution V = V[w] of the equation 
— A x V[w] = n[w], x € IR 3 , satisfies 

II WH lUepRS) < C||nH|| L2(R3) . (2.4) 

Proof. Since V = — ■— i * n, we have W = . f n * n, and the estimate follows from the 
generalized Young inequality. □ 

Remark 2.4 Note that n G L 2 (IR 3 ) does not yield (via the generalized Young inequality) 
a control of V in any U — space (even n G L P (IR 3 ) with p G (3/2,2] would not "help"). 
However, the operator Q[V] involves only the function 6V, which is slightly "better behaved". 
We anticipate that we will later recover some information on the potential V via new a-priori 
estimates on the electric field W[w] (see Corollary \4- IT}) . 

Omitting the time-dependence we have 

5V{x,r,) = V(x + -)-V(x- -) = —j ^ ^ d£ 

f f(y,v) n l w ]( x - y)dy, 



4-7T 
R 3 



with the "dipole-kernel" f(y; n) := 
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Proposition 2.5 For all w G X and fixed rj G IR 3 , we have 6V[w](., 77) G L 9 (IR 3 ), 6 < g < 
oo. Moreover 

\\SV[w](.,r,)\\ LOO{Ml) < Cl^l^llnHIL^Ra)- (2-5) 
Proof. By using the triangle inequality, 

M 



|y-3l|y + 3l' 



and the transformation y = \n\x, we estimate for 3/2 < p < 3 



where e G K 3 is some unit vector (due to the rotational symmetry of ||/(.; ^II^phr 3 ) with 
respect to rj). Young inequality then gives SV(.,rj) G L <? (1R 3 ), 6 < q < oo, and the assertion 
holds. □ 

In most of the literature the Wigner operator is defined on L 2 (IR^) for bounded potentials 
V, cf. |Mr1IME1IA7TD] . For our nonlinear problem l[T"T |) - (fi3|L however, V G L°°(1R 3 ) does 
not hold. As a compensation we shall hence exploit the additional regularity of the Wigner 
function to define the quadratic term G[V[io]]tu (cf. Prop. 2.8 in |Maj for a similar strategy). 

Proposition 2.6 Let u G X and V v u £ X be given. Then, the linear operator 

z\ — >0[V[z]]u, 

with the function V[z] = — * n[z], is bounded from the space X into itself and satisfies 

\\e[V[z]]u\\ x < C{\\u\\ x + \\V v u\\x}\\z\\x, Vz G X. (2.6) 

Proof. To estimate ||0[V[.z]]u||x we shall consider separately the two terms of the equivalent 
norm 

3 

2 |i.i|2 , ||„.2„.||2 



\U\\ X := \\U\\2 



+ £II«MI2- ( 2 - 7 ) 



i=l 

First, by denoting u := JF^-^u, we get 



|e[f[*]]u||! = JJ\8(V[z])(x,n)u(x,n)\ 2 dxdr, < Jj\\5(V[z]){.,v)\\Uu(x,n)\ 2 drjdx 
< C\\z\\ x [[(\r,\V 2 \u(x,r,)\y drjdx < C\\zf x (\\u\\ 2 + \\V v u\\ 2 ) , (2i 



by applying first the Plancherel Theorem, then Holder's inequality in the x variable, the 
estimates (|2.3j) . (|2.5(l for the function <5V[z], and finally, Young inequality and the Plancherel 
Theorem to the last integral. 
For the second term of ||©[V[,z]]ii||^ we shall use 

v 2 Q[V]w(x,v) = -e[^V]w(x,v) + n[diV}(v iW )(x,v) + @[V]v 2 w(x,v), (2.9) 

with the pseudo-differential operator 

Q[V] := i(6 + V)(x,^y (S + V)(x, V ) := v(x + |) + v(x - fj . (2.10) 
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Here and in the sequel we use the abreviation di := d Xi . Q2.9JI is now estimated: 

\\vfe[V[z]]u\\ 2 < i||<5(^y[z])n|| 2 + ||5 + (a i y[z])^n|| 2 + || ( 5y[z]a^|| 2 (2.11) 
The first two terms of (|2.11j) can be estimated as follows: 

II<5(^MHIl2(ir 6 ) < 2 ll^ 1/ WllL2(R3 ) ||n|| L2(]R 3. L « )(]R 3 )) 
< CHzllxIKl + H^ulliafRe), 



by applying Holder's inequality, (j2.3j) and the Sobolev imbedding u(x,.) G -ff 2 (IRf ? ) <— > 

l°°0r3). 

< C||z||x||(l + ^ 2 )n|| 2 , (2.12) 

by the Sobolev imbedding and V v u(x, . ) G iJ 1 ^^) ^ L 4 (IR^), and by estimate (|2~i|) for 
W[z] and (|2~3|). For the last term of (|2~TT|) we estimate as in (|2~H)l : 



< C||z|ft [f(\ V \V 2 %u(x,T,)) 2 dr,dx 



< c\\z\&(\\%ni+\\v%ni) 

< c\\zf x (wdim +11^(^)111+11^^11) 

< C\\zf x (\\(l + vf)u\\l + \\vfV v u\\l), 



by interpolation and integration by parts. 

This concludes the proof of estimate (|2.6|) . rj 
Remark 2.7 T/ie previous proposition shows that the bilinear map 

(z,u) i — ► 8[V[js]]ti 

is well-defined for all z,u G X, subject to V v u G X. The unusual feature of the above 
proposition is the boundedness of this map with respect to the function z appearing in the 
self- consistent potential V[z]. This is in contrast to most of the existing literature (I A (JE\ 
\MBl \Mldj ). where the boundedness of the pseudo-differential operator Q\V[z]] (with z fixed) 
is used. However, this can only hold for bounded potentials V . 

2.2 Dissipativity of the linear equation 

In our subsequent analysis we shall first consider the linear Wigner-Fokker-Planck equation, 
i.e. equation Ql.l|) with V = 0. The generator of this evolution problem is the unbounded 
linear operator A : D{A) — ► X, 

Au := — v ■ V x u + (3div v (vu) + aA v u + 2^div v (V x u) + aA x u, (2-13) 

defined on 

D(A) = {u G X | v ■ V x u,v ■ V v u, A v u,div v V x u, A x u G X}. ( 2 -14) 
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Clearly, Cq°(IR 6 ) C D(A). Hence, D(A) is dense in X. Next we study whether the operator 
A is dissipative on the (real) Hilbert space X, i.e. if 

<Au,u>x< 0, VtieD(i) (2.15) 

holds. 

Lemma 2.8 Let the coefficients of the operator A satisfy ao~ > 7 2 . Then A — kI with 

k := -13 + 9a (2.16) 

is dissipative in X . 

The proof is lengthy but straigthforward and deferred to the Appendix. 



By Theorem 1.4.5b of |Paj its closure, A — kI = A — kI is also dissipative. 
A straightforward calculation using integrations by parts yields 

<Au,w>x = <u,A\w>x + <u,A2W>x, \/u,w&D(A), 

with 

A\w = v • V x w — j3v ■ V v w + aA v w + 2^div v (V x w) + aA x w, 
<u,A* 2 w> ji = ^(-^(3jJvtwu + ^aJJvfw Vz u + YO-ffvfwu+^jffvfw Xz u 



i=l 



Hence, A* \d{a) ~ the restriction of the adjoint of the operator A to D(A) - is given by 
A* = A\ + A\. A* is densly defined on D(A*) 5 D(A), and hence A is a closable operator 
(cf. Theorem VIILl.b of jRHU). Its closure A satisfies (A)* = A* (cf. |Ri5T| . Theorem 
VIII.l.c). 

Since < A*u, u > = < Au, u > the following lemma on the dissipativity of the operator A* 
restricted to D(A) holds. 

Lemma 2.9 Let the coefficients of the operator A satisfy aa > r y 2 . Then A*\ D ^ —kI is 
dissipative (with k as in \2.16\) ). 

Next we consider the dissipativity of this operator on its proper domain D(A*), which, 
however, is not known explicitly. To this end we shall use the following technical lemma, 
which, for a matter of generality, is stated in the space with the symmetric x,t>-weight. The 
proof is defered to the appendix: the arguments are inspired by |ACD| . |ASj . but there are 
also similar results for FP-type operators in |HelNl iHerNj . e.g. 

Lemma 2.10 Let P = p(x, v, V x ., V^) where p is a quadratic polynomial and 

D{P) := C °°(1R 6 )CZ := L 2 (1R 6 ; (1 + \x\ 2 + \v\ 2 ) 2 dx dv). 

Then P is the maximum extension of P in the sense that 

D(P) := {u G Z | the distribution Pu 6 Z} . 
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We now apply Lemma l2.1Ul to P = A* — kI, which is dissipative on D{P) C D{A). Since A* 
is closed, we have D(A*) = D(P) = {u G X \ A*u G X} and A* — kI is dissipative on all of 
D(A*). 

Applying Corollary 1.4.4 of |Paj to A—nl (with (A)* = A*), then implies that A—kI generates 
a Co semigroup of contractions on X, and the Co semigroup generated by A satisfies 

||e i3 n||^ < e Kt |H|^, ueX, i > 0. 
Since \\.\\x and \\-\\ x are equivalent norms in X with 

\\u\\ x < \\u\\ x < 4||u||_y, 



we have 



e tA u\\ x < 4:e Kt \\u\\ x , u£X, t>0. (2.17) 



3 Existence of the local- in- time solution 

In this section we shall use a contractive fixed point map to establish a local solution of the 
WPFP system. To this end the parabolic regularization of the linear WFP equation will be 
crucial to define the self-consistent potential term. 

3.1 The linear equation 

First let us consider the linear equation 

w t = Aw(t), t>0, w(t = 0) = wq G X. (3.1) 
By the discussion in Subsection 12.21 its unique solution w(t) = e wq satisfies 

IK*)lk < 4e Kt |KHx, Vt>0. (3.2) 
Actually, the solution of the equation can be expressed as 

w(x,v,t) = wq(xq, vo)G(t, x, v, xq, vq) dxQ dvo, V(x,u)GlR 6 , (3-3) 



where the Green's function G satisfies (in a weak sense) the equation (|3.1(l and the initial 
condition 

lim G(t, X, V, Xq, Vq) = 5(x — Xo , v — V o) , 

for any fixed (a; ,^o) G 1R 6 (cf. Def. 2.1 and Prop. 3.1 in |SCI)M| ). 
The Green's function reads 

G(t,x,v,x ,v ) = e m F(t,X_ t (x,v) -x ,X_ t (x,v) -v ), (3.4) 

with 

F(txv) - I f v(t)\x\z + \(t)\vf + ^t)(x.v) \ 

K ' ' ' (2vr)3(4A(tMt)-^ 2 W) 3/2 I m>{t)-H 2 {t) ]' 

The characteristic flow $>t(x,v) = [Xt(x,v), Xt(x,v)] of the first order part of (|2.13|) . is given 
for (3 > by 



X t (x,v) = x + v 
X t (x,v) = ve-f*, 



9 



and $t{x, v ) = [x+vt, v], for (3 = 0. The asymptotic behaviour of the functions X(t), v{t), fi(t) 
for small t is described (also for (3 = 0) by 



X(t) 
u{t) 

A*(t) 
And hence: 



at + a 

e 2/3t_! 



e 2,3t„ 4e fit + 3 j 
2^ h ^ 1 



+ 7 



2« 







2 ( e m 



i) 



a 



a 



2/3 
(^) 

/(*) 



2 . 2(l-e' 3t ) 



at, t 
at, t 

-27*, « 



0. 
0, 
0. 



4A(t)i/(t) -// 2 (*) 



4(aa - 7 2 )t 2 > 



With these preliminaries, the following parabolic reguralization result can be deduced. 

Proposition 3.1 For each parameter set {a, (3, 7,cr}, there exist two constants 

B = B(a, (3, 7, a) and Tq = Tq(ol, /3, 7, a), such that the solution of the linear equation h'J.l\) 

satisfies 



\V v w 

\V T W 



< st-vv'iKlk, 



V0 < t < T , 
V0 < t < T , 



(3.5) 
(3.6) 



for all wo G X. 

The proof is similar to |Carj and it will be defered to the Appendix. 



Remark 3.2 (a) Observe that the functions V^w, V v w G C((0, oo);X). The local bounded- 
ness of V x w, V v w on any interval (r, r + To] follows from \3. ty) and Prop. \3.1\ 
(b) Note that the strategy of the next section will not work in the degenerated parabolic case 
aa — 7 2 = 0, since the decay rates of Prop. I.V.il would then be t~ 3 / 2 , which is not integrable 
at t = 0. Alternative strategies for this degenerate case were studied in I ALMS)! . 



3.2 The non-linear equation: local solution 

Our aim is to solve the following non-linear initial value problem 

w t (t) = Aw(t) + Q[V[w(t)]]w(t), Vt>0, w(t = 0) = w £X, (3.7) 

where the pseudo-differential operator O is formally defined by (|1.2j) and the potential y[io(t)] 
is the (Newton potential) solution of the Poisson equation 

-A x V{t,x) = n[w(t)](x) = I w{t,x,v)dv, x G IR 3 , (3.8) 

JlR 3 

for all t > 0. Actually, if we assume w(t) G X for all t > 0, then the function n[w;(t)] is 
well-defined for all t > (cf. Prop. 12.1(1 . and the solution V[io(t)] satisfies the estimates of 
Propositions for all t > 0. 

The Propositions 12.61 and 13.11 motivate the definition of the Banach space 

Y T := I z G C([0, T]; X) \ V v z G C((0, T];X) with \\V v z(t)\\ x < Ct^ 2 for t G (0, T)}, 

endowed with the norm 

||z||y T := sup ||z(t)||x+ sup \\t 1/2 V v z(t)\\ x , 
te[o,T] te[o,T] 
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for every fixed < T < 00. We shall obtain the (local- in-time) well-posedness result for the 
problem 1)3.7)) by introducing a non- linear iteration in the space Yt, with an appropriate 
(small enough) T. 

For a given w G Yt we shall now consider the linear Cauchy problem for the function z, 

z t = Az(t) + e[V[z(t)]]w(t), VtE(0,T], z(t = 0) = w £X, (3.9) 

with < T < Tq and To is defined in Prop. 13.11 According to Prop. 12. 61 the (time-dependent) 
operator Q[V [. }}w(t) is, for each fixed t G (0, To], a well-defined, linear and bounded map on 
X, which we shall consider as a perturbation of the operator A. 

Lemma 3.3 For all wq G X and w G Yt, with T < Tq, the initial value problem 

z t = Az(t) + e[V[z(t)]]w(t), VtG(0,T], z(t = 0) = w , 
has a unique mild solution z G C([0, T];X), which satisfies 

z(t) = e a w + f t e {t ~ s) ' J e[V[z(s)}}w(s)ds, Wg[0,T]. (3.10) 
J 

Moreover, the solution z belongs to the space Yt- 

Proof. The first assertion follows directly by applying (a trivial extension of) Thm. 6.1.2 
in jPij: 

For any fixed w G Yt, the function g{t, .) := G[y[.]]^(t) is a bounded linear operator on 
X for all t G (0,T), and it satisfies g G T 1 ((0, T);B{X)) n C((0, T];B{X)) (by Prop. l2~6l) . 
Moreover, by estimates ()2.17)) . ()2.6)) . the following inequalities hold 



\\z(t)\\x < 4e K *||^o||x + 4 / e< t -^C{\\w{s)\\x + \\yvw{8)\\x}\\z{8)\\ x ds (3.11) 

Jo 

< te Kt \\ W(s \\x + ±Ce KT \\w\\Y T I {l + s^'^Wz^Wxds, (3.12) 

J 

for all i G [0, T]. Then, by Gronwall's Lemma, 

\\z{t)\\ x < 4e- T || Wo |U [l + 4C||^||y- T e( K ' r+4C7eKr l^ll^^ +2:rl/2 ))(t + 2^/2)] ^ (3J3) 
for all t G [0, T\. By differentiating equation ()3.1Uj) in the f-direction, we obtain 

V v z(t) = V v e tA w + V v e {t - s)A g(s,z(s))ds, VtG[0,T]. (3.14) 

J 

Using the estimates 1)3.5)1 . 1)2.6)1 . and ()3.13)) then yields 

||V^(*)|U < Bt-vv^i^olk 

+ B\\w\\ Yt [\t - s)" 1/2 e K( '- s) C{ 1 + s^ 1/2 } ||z(s)||x rfs 
Jo 



< Br 1/2 e Kt \\w \\x + 45Ce 2KT ||u;o||x||u;||y T 



vr + 2t 



1/2 



+ 4C|| W ||y T e( KT+4CeKT ll u 'll^(^ 1/2 )) Ut 1 / 2 + -vrt + -t 3 / 2 

V 2 3 



(3.15) 



for all t G [0, T]. The continuity in time of V v z can be derived from ()3.14)) by using Remark 
13. 21 and the fact that g(t, z{t)) G C((0, T];X). Hence, the function z belongs to the space Yt- □ 
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We now define the (affine) linear map M on Yj> (for any fixed < T < To): 

Mw : = z. 



w 



where z is the unique mild solution of the initial value problem 1)3.9(1 . According to Lemma 
13.31 z eYt- Next we shall show that M is a strict contraction on a closed subset of Yt, for T 
sufficiently small. This will yield the local- in-time solution of the non- linear equation (|3.7|) . 

Proposition 3.4 For any fixed wq G X, let R > m&x{4:, B}e K \\wo\\x , with the constant B 
defined in Prop. \3.1\ Then there exists a r := t(||wo||x, B) > such that the map M, 

— rt — 

(Mw){t) = e tA w + e {t - s)A @[V[Mw{s)]]w(s)ds, Vt€[0,r], (3.16) 
Jo 

is a strict contraction from the ball of radius RofY T into itself. 

Proof. By (the proof of) Lemma 1531 the function z = Mw G Y T satisfies (|3.13|) . Under the 
assumption ||«j||y T < R, this estimate reads 



\\Mw(t)\\ x < 4e KT |K||x 
If we assume 



1 + 4Ci?e( KT+4CReKT ( T+2rl/2 ))(t + 2i 1 /2)l , Vt G [0, r]. 



4e KT \\w \\ x 



1 + 4CRe( KT+ACReKT ^ +2Tl/ ^)(r + 2T 1 ' 2 ) 



< 



R 



(3.17) 



then \\Mw{t)\\ x < f • Similar to (ETXKl) we have 

\\V v Mw(t)\\ x < Bt^ 2 e Kt \\w \\ x + 4BCRe 2KT \ 



w \\x 



+ 4CRe( KT+ACReKT ^ +2Tl/ ^) ( At 1 ' 2 + ±7Ct + |t 3/2 



If we assume 

Be KT \\w \\x + 4BCRe 2KT \\w \\x 



vrr 1/2 + 2r + 



+ 4C i2e (-+4Cite-(r+2r 1 / 2 )) ( 4t + ^t 



4 
3 



< 



i? 



(3.18) 



then 



R 



t l l 2 \\V v Mw(t)\\ x < -, Vt G [0,r] 



Let us now choose 



min < 1, 



/ i?/3-4e K |K|| x 
\48CR\\w \\ x e 2K + 12CReK J ' 

R/3 - Be K \\w \\ x 



4BCRe 2K \\w \\ x vr + 2 + ACRe( K + 12CReK )(^Tr + 



(3.19) 



which is positive since max{4, i?}e K ||i(Jo||x < R- Then, the estimates (|3.17() and (|3.18|) hold, 
and hence the operator M maps the ball of radius R of Y T into itself. 
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To prove contractivity we shall estimate ||Mit — Mw;||y T for all u, w £ Y T with ||«||y T , IMIvv < 
R. Since 



Mu(t) - Mw(t) 



^- s ) A Q[V[(Mu - Mw){s)])u{s) ds 



+ / e {t - s)A e[V[Mw(s)}]{u-w)(s)ds, Vt G [0,r], 
Jo 



by analogous estimates, 

\\Mu{t) - Mw{t)\\ x < 4CRe KT { \ "(1 + s" 1/2 ) || (Mu - Mw)(s) \\ x ds 



+ \\U - W\\Y, 

and, by applying Gronwall's Lemma: 



\\Mu(t) - Mw(t)\\ x < 4Ci?e* 



t + 2t 1/2 + 



+ 4C7 J Re( Kr+4CfleKT ( T+2rl/2 )) [ 2t + 2t 3 / 2 + V 



u - w\\y T , Vi G [0,r]. 

(3.20) 



t 1 / 2 ^ — w\\y t ■ 



By using < t < r < 1, we obtain 

||M«(f) - Mw(f)||x < 4Ci?e K [3 + 18C7iie (K+12Cfle,t) 
Similarly, 

\\V v Mu{t) - V v Mw(t)\\ x < CBRe kT l^j\t - s)'^ 2 {l + s" 1 / 2 ) ||(M« - Mw)(s)\\ x ds 

+ At-sJ-Va^ + a-i/a)^ 

and, by using estimate (|3.2()|) . 



1 + 4CRe K ( 3 + 18C7?e( K+12Cife ^ r 1 / 2 



l + 4Ci?e M 



\\V v Mu(t) -V v Mw(t)\\ x < CBRe 

• (vr + 2t 1 / 2 ) \\u-w\\ Yt , Vt€[0,r]. 

Then, by exploiting < r < 1, 

t 1/2 \\V v Mu(t) - V v Mw(t)\\ x < CBRe K (n + 2 

• (3 + l8CRe^ +12CRe ^ 
When choosing r > small enough, estimates ()3.2Uj) . H3.21j) imply 

\\Mu - Mw\\cQp, T ] ; x) < C\\u - w\\ C ([o, T ];X), 

for some C < 1, and the assertion is proved. rj 

Corollary 3.5 There exists a i max < 00 such that the initial value problem {3. 1\) has a 
unique mild solution w in Yt, VT < i max , which satisfies 



t 1/2 \\u- w\\ Yt . (3.21) 



w (t) = e tA w + e^ A @[V[w(s)}]w(s)ds, VtG[0,T]. 
Jo 



(3.22) 



Moreover, i/t max < 00, then 



lim ||u;(i)||x = 00. 
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Proof. The solution of the problem is the fixed point of the map M previously introduced. 

By Prop. 13.41 this solution exists for a time interval of length r (depending only on ||wo||jf) 

and it belongs to the space Y T . Since, in particular, w{t) G X, the solution can be repeatedly 

continued up to the maximal time t max . It will then belong to Yt, VT < i max . 

If the second assertion of the corollary would not hold, there would be a sequence of times 

tn T t max such that ||w(t n )||x — C for all n. Then, by solving a problem with the initial value 

w(t n ), with t n sufficiently close to t max , we would extend the solution up to a certain time 

t n + T(\\w(t n )\\x) > imax- This construction would contradict our definition of i max . 

The uniqueness of the mild solution follows by arguments analogous to those in the proof of 

Thm. 6.1.4 in jPa]. n 

Remark 3.6 Note that the last statement in the thesis of the Corollarv VJ. 51 differs from the 
standard setting (cf. Thm. 6.1. 4 in \PcHj ). For i max < oo we conclude the 'explosion' of 
max i n X and not only in Yt- This is due to the parabolic regularization of the 
problem \3.T\j . 



4 Global-in-time solution, a-priori estimates 

In this section we shall exploit dispersive effects of the free transport equation to derive an 
a-priori estimate on the electric field. This is the key ingredient for proving the main result 
of the paper, the global solution for the WPFP system: 

Theorem 4.1 Let Wq G X satisfy for some lo G [0, 1) 



woix — $(t)v, v) dv 



£6/5 (]R 3) 



< Ct#(*)~ 



ViG(0,T], VT>0, 



with #(t) := ±=|—- for f3 > 0, and 0(t) = t for (3 = 0. Then the WPFP equation admits 
a unique global-in-time mild solution w G Yt, V0 < T < oo. 

In order to prove that t max = oo , we have to show that ||u;(i)||x is finite for all t > (cf. 
Corollary 13 .5j) . To this end, we shall derive a-priori estimates for ||w;(i)||2 and | 2 w(t)||2. 
Thus, the proof of Thm. 14.11 will be a consequence of a series of Lemmata, in particular 
of Lemma 14.21 and Lemma 14.191 In the sequel, w(t) denotes the unique mild solution for 
< t < T, for an arbitrary < T < t max . 

Lemma 4.2 For all wq G X , the mild solution of the WPFP equation \3.T\j satisfies 

\\w(t)g < e 3/3t |Kll2, yt G [0,T]. (4.1) 

Proof. Roughly speaking, this follows from the dissipativity of the operator A — ^§-1 in 
L 2 (IR 6 ) (cf. (|5.2|) ) and the skew-symmetry of the pseudo-differential operator. However, 
since we are dealing only with the mild solution of the equation, the proof requires an 
approximation of w by classical solutions. 

Since the solution satisfies w G Yt, VT < i max , Prop. 12.61 shows that the function f(t) := 
e[V[w(t)]]w(t), t G (0,t max ) is well defined and it is in C((0, t max ); X) n ^((0, T); X), V0 < 

T < imax' 

For < T < i max fixed, let us consider the following linear inhomogeneous problem: 

^-y(t) = Ay(t) + f(t), te[0,T], y(t = 0) = wq G X. (4.2) 
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Its mild solution in [0,T] is the function w, due to the uniqueness of the mild solution of 
problem (|3.7jh For this linear problem, we can apply Thm. 4.2.7 of |Paj : The mild solution w 
is the uniform limit (on [0, T]) of classical solutions of problem (|4.'2[) . More precisely, there is a 
sequence {^olneiN C £>(^4), w>g — ► wg i n -X", and a sequence {/n(*)} C C 1 ([0, T];X), f n (t) — > 
/(i) in L 1 ((0, T);X). And the classical solutions y n G C 1 ([0,T];X) of the corresponding 
problems 

j t y n (t) = Ay n (t) + f n (t), te[0,T], y n {t = 0) = teg, (4.3) 
converge in C([0,T];X) to the solution «; of problem (|4.2j) . 

We shall need these approximating classical solutions y n in order to justify the derivation of 
the a-priori estimate: Multiplying both sides of (|4..3[) by y n (t) and integrating yields 

^l|yn(t)||! < Y\\yn{t)\\l + JJyn{t)fn{t)dxdv, 

since the operator A — ^ is dissipative in L 2 (IR 6 ) (cf. (j5.2|> ). By integrating in i and letting 
n — ► co, we have 

IK*)lli < \\wo\\l + 3pJjw(s)\\lds + 2j^JJw(s)f(s)dxdvds, VtG[0,T]. 

The second integral is equal to zero since the pseudo-differential operator O is skew-symmetric. 
Hence, applying GronwalFs Lemma yields 

\\w(t)\\l < e^l^olli, VtG[0,T]. (4.4) 

□ 

In order to recover similar estimates for || \v | 2 u>(t)||2, we first need a-priori bounds for the 
self-consistent field E = VV. To this end, we are going to exploit dispersive effects of 
the free streaming operator. We shall adapt to the Wigner-Poisson and Wigner-Poisson- 
Fokker-Planck problems the strategies introduced for the (classical) Vlasov-Poisson problem 
f lLTlIPe] ). and for the Vlasov-Poisson-Fokker-Planck problem f fBoll lBo2l |Ha2] 1. 

4.1 A-priori estimates for the electric field: the Wigner-Poisson 

case 

To explain the strategy, we first consider the (simpler) WP problem: Let us assume that w wp 
is a "regular" solution of the WP problem (e.g., let w vp (t) £ L 2 x (H^),V x V[w vp }(t) £ Cb(M 3 ), 
uniformly on bounded t-intervals) for which the Duhamel formula holds: 

w wp (x,u,t) = w^(x-tv,v) + j (eiVlw^w^Yx - sv,v,t- s)ds. 

We formally integrate in v: 

n[w; wp ](x,t) = / Wn p {x -tv,v)dv + t ( ( 9[y[w; wp ]]u; wp ) (x - sv, v, t - s) dv ds 
Jm 3 Jo Jn* v J 

=: ri?(x,t) + r£?{x,t), 
and split the self-consistent field accordingly: 

E- p (x,t) := \^L* x nQ P (x,t) = * x J w?(x - tv, v) dv, (4.5) 

E™ p (x,t) := * x J^J(e[V[w vp ]}ur P Yx-sv,v,t-s)dvds, (4.6) 
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with A = ^. 

Then, we can estimate separately the two terms E^ v {t), E™ p (t) by exploting the properties 
of the convolution kernel l/\x\ (cf. |LPl EHJ for VP, [BoTl IEo21 IHa2j for VPFP, |ALMS| for 
WPFP). To this end, we need an appropriate redefinition of the pseudo-differential operator 
G[V] in p.2|) . It is inspired by the operator V x V-V v w in the VP equation that can be 
recovered from 0[V]u; in the semiclassical limit (cf. Remark 14.5(1 . 

Let us recall that Q[V]w(x,v) = T~\ v (i 5V(x, r])^F v ^ ri w(x, r])j . We can rewrite 

x+n/2 1/2 
5V(x,rj) = J V x V(z)-dz = J r] ■ V X V (x - rrj) dr = rj-W(x,r]), (4.7) 
x-77/2 -1/2 

with the vector-valued function 

1/2 

W(x,rj) := J V x V(x - rr])dr , V (x, rj) G M 6 . 
-1/2 

Then, we define the vector- valued operator 

r[V x V]u(x,v) := F-± v (w{x,r))F v ^ n u{x,T])). (4.8) 

It holds: 

Lemma 4.3 Let V X V £ C B {M 3 ). Then 

1. W(x, V ) GC B (IR 6 ), II^IU < ||V x ^||oo; 

2. T[V X V] : L 2 (IR 6 ) -> L 2 (IR 6 ) and, for all u G L 2 (IR 6 ) , 

HrfV^^uH^^e) < ||V a .V r ||oo||«|li2(]Re); 

3. F[V X V] : L 2 x (Hl) L 2 x (Hl) and, for all u G L 2 (i?i), 

H r [VxVlu||i5(Hi) < ||V a .V r || QO ||«||^( fl i). (4.9) 
Proof. The first and the second assertion are obvious. For (|4,9|) we use 

d v T[V x V]u(x,v) = ^^(rijWix,!])^ ^rjU^x,!])^ = T[V x V]d Vj u; j = 1,2,3. (4.10) 

□ 

Lemma 4.4 Let V X V G Cb(IR 3 ) and u G L 2 x (Hl). Then 

@[V]u(x,v) = dw v (T[V x V]u)(x,v) (4.11) 
Proof. By the definition (|4.7j) and Lemma [Ql 

ll^c^lloo < MII^C^Iloo < Nl|v x y||oo. 

Thus, ||@[V r ]'u|| i 2(] R 6s < || Va; V||oo ||it||x,2 (^-ij; the right hand side of equation (|4,11[) is also 
well-defined in L 2 (IR 6 ) by estimate ([4.9)1 . Equality then follows by equation (|4.7|) and 

3 

iF-i v (v-W(x,r))F v ^ v u{x,ri)) = ^^jl^xVHx^) = dir v (T[V x V]u) (x,v). 

□ 
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Remark 4.5 (The semiclassical limit) The correctly scaled version of the pseudo-diffe- 
rential operator with the reduced Planck constant fi = reads 



O h [V]w(x,v) 



(2tt) 3 / 2 



V{x + §7?) - V(x - §77) 



n 



Under the assumptions of Lemma \4-4[ we thus have 



F v ^ v (Q n [V]w(x,v)) = -5V{x,hn)J r v ^ rj w(x,n) 

= iW(x,hn) • r]F v - ¥ r)W(x,7)). 

The limit fi — > then yields: 

iW(x,hn) • nF v ^ v w(x,rf) — > iV x V(x) ■ r) F v ^ n w{x, rj) = F~}> v (y x V(x) ■ V v w(x,v)j; 
and hence 

& h [V]w(x,v) — ► V x V{x) -V v w{x,v) m L 2 (E 6 ), 
which is the non-linear term in the VP equation. 

Using the redefinition l|4.11J) of the pseudo-differential operator, and under the additional 
assumptions w wp £ H*(L%), AF[w wp ] £ C B (1R 3 ), we have for s S E 

(e[U[tt; wp ]]u; wp )(2; - sv,v) = div^(V[V x y[w wp ]]u; wp (x - sv, vfj 

+ s div x (t[V x V[w^}]w^ (x - sv, v). (4.12) 

Thus, the field EJ P in Q4.6|l can be rewritten as (j = 1, 2, 3) 



{E?)j(x,t) := A|— |g * x div x J sj (TjV x y[u/ wp ]]™ wp J (x - sv,v,t - s)dvds 

= A ]T l SjklX{2 ** f* J (r k \y x V[w~>]]vr») (x - sv, v,t-s) dv ds. 

k=l ^ 

The following two lemmata are concerned with giving a meaning to the definition (|4.13l) of 
the field E±, independently of the previous derivation. 

Lemma 4.6 For all u £ L 2 (IR 6 ) and E £ L 2 (IR 3 ) the following estimate holds 

< Cs~ 3/2 ||S|| 2 ||n||2, Vs>0. (4.14) 



(T[E]u) (x — sv, v) dv 



Proof. Since the operator T[. \ was originally defined for E £ Cb(IR 3 ), we shall first derive 
(j4.14j) for E £ Cq°(IR 3 ) and conclude by a density argument. 

By the definition (|4.8j) and by several changes of variables, the following chain of equalities 
holds: 

(T[E]u)(x,v) = (2tt) 3 / 2 [r£v(W(x,rj))* v u] (x,v) 

1/2 1/2 

E{x — rrf) e tr] ' z dr dnu(x,v — z) dz = // / y-r^E{x — rf) e tf, '~ dr dfju{x,v — z) dz 



-1/2 



-1/2 



J J E(x — fj) e lf, ' z dfj J u(x, v — rz) dr dz 

-1/2 
1/2 

(2vr) 3 / 2 j^- z E{z) J u(x,v-rz) 



1/2 

E(-fj) e tfl ' z dfj e tx ' z J u(x, v - rz) dr dz 

-1/2 



dr e lx ' z dz. 



-1/2 
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Hence 



J (T[E]u) (x - sv, v) dv = (2tt)§ J T^~ z E{z) ( j J*\ 



(2vrs) E 



u(x - sv,v - rz) dr e lsv ' z dv \e lx ' z dz 



V2 v \ . _ 

u(x - v , rz) dr e tx ' z dz . 

-1/2 S 



Then, 
J(T[E]u) (x - sv,v) dv 



L*(M 3 X ) 



< 



< 



\\Eh 


rl/2 

/ T x 

J -1/2 


(2vrs) 3 


\\E\\ 2 


in 


(2vrs) 3 


17-1/2 



1/2 / f \ 

T v ->z u(x — v , rz) 

^1/2 V S J 



L2(R|xKi.) 



by applying Holder's inequality first in the z integral and then in the r integral. Finally, it 
remains to prove that 



1/2 



\\T v ^ z (u{x -v, rz)) | 

1/2 v s ' 



3„Tr>3\dr 



L 2 (fRlxTR 3 2 ) 



s 3 \\u(x,v)\\ 2 L2{nlx]R3vy 



This is obtained by using repeatedly Plancherel's equality: 



1/2 



1/2 



u(x — v, rz) 



-1/2 



1/2 

/' 

-1/2 

s 3 \\u(x,v)\\ 2 L 2 



J \\T x ^(s z e- is ^ rz T v . 



h dr = 

-1/2 
, s ^ +z )U{x,V 



e iv ^u(x, rz 

v s 



Si* 



1/2 



Si* 



dr 



dr 



s,* 



-1/2 



□ 

Remark 4.7 Observe that the exponent of the variable s recovered in the Lemma is the 
same as obtained for the VP case (cf. }Pej /. e.g.) in the I? -estimate of Eu(x — sv, v) dv. 
In the classical case, analogous LP -estimates hold in addition. In the quantum counterpart, 
instead, the L 2 -framework is the only possible, since the estimate had to be derived in Fourier 
space. Moreover, to derive a more refined version of this basic estimate (cf. \LP\ \BoT^ ). the 
non-negativity of the classical distribution is a crucial ingredient. And this non-negativity 
does not hold for Wigner functions. 

The following lemma is an immediate consequence of Lemma [4.61 We shall need the notation 
Vt*, ■= {E G C((0,T]; Ll(JR 3 ) \ \\E\\ Vt ^ < oo} 

with 



\E\Wn 



T,u 



sup r\\E(t)\\ L 2. 

0<t<T 



Lemma 4.8 For any fixed T > 0, let w G C([0, T]; L 2 v ), and let wq satisfy for some to G 
[0,1): 



wq(x — tv, v) dv 



< C T r", Vt G (0,T]. 



(B) 
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Then, there exists a unique vector-field E £ V T 1 which satisfies the linear equation 



3 o„_. , x.. UI2 



= a — 3 us ** y o s y (rfc[Eo + EVj(x-s^u,t- s )^ds; i = 1,2,3 

(4.16) 



fc=l 

with Eq defined by (A = ^ 

E(j(x,t) := A, — pr *x / wo(x — tv , v) dv 



Proof. Q4.16JI has the structure of a Volterra integral equation of the second kind. Hence, 
we define the (affine) map M : V TuJ _i — > V Tui _i by 



_> 

(ME) j (x, t) := A J2 ~ 3XlXk ^ 5 6 ^ 2 * x J* 8 J (r k [E + E]w) (x - sv, v,t - s) dv ds. 
Applying the generalized Young inequality to the definition of Eq yields 



\\E (t)\\ 2 < C\\J w (x-vt,v)dv\\ L6/5{nl) , Vie(0,T]. (4.17) 
Thus, by Lemma 14.61 the second convolution factor in (|4.16|) is well-defined and 

[ (r k [E +E]w)(x-sv,v,t-s)dv r2 ,. < Cs- 3 / 2 \\(E +E)(t-s)\\ 2 \\w(t-s)\\ 2 ,Vse (0,t]. 

By classical properties of the convolution with m (cf . [S3 ) and the Young inequality, we get 

||(Mi%(t)|| 2 < C fl-(\\E G {t-s)\\ 2 + \\E{t-s)\\ 2 )\\w{t-s)\\ 2 ds, Vie(0,T]. (4.18) 
Jo V s 

Hence, the map M is well-defined from V Tui _i into itself and satisfies 

||ATE(t)|| 2 < C(C T + sup s"-3\\E{s)\\ 2 ) sup \\w(s)\\ 2 (t 1 ^ + , Vie(0,T]. 

v se(o,T] y se[o,T] v y 

Since the map is affine, we have (by induction) for all t E (0, T] 

||M n £(i)-M n £(t)|| 2 < C sup ||te(a)|| 2 f \\M n - x E(s) - M n ~ l E(s)\\ 2 ds 

se[o,T] Jo Vt- s 

< (C sup Hs)|| 2 ) C n _i / ^5= ds SU P (^ll^-^OOIk), 



with 



n+i , , „ / 1 n + 2 
ds = t~ — "B -, w 



\2' 2 



where B denotes the Beta function and T the Gamma function. Thus, the map M n is con- 
tractive for n large enough and admits a unique fixed point E E V T w _ i . rj 

With = the above lemma yields the regularity of the self-consistent field in the 
WP equation: It satisfies V x y[w wp ] = E™ p + Eq p E V Tu) _i, under the assumptions that 

w wp E C([0,T];L 2 xv ) and Wq p satisfies (Bj). 
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Proposition 4.9 For any fixed T > 0, let w wp 6 C([0, T]; L 2 , v ) be a mild solution of the WP 



equation with \\w wp 



Wn 



I2, and with the initial value Wq p satisfying condition iTb)) , 



Then, the self- consistent field satisfies the following estimates for all t £ (0, T] : 



< C 



Wq p (x — vt, v) dv 



< cc T t-", 



Vll2>SUp se ( ,T] i S 



Wq p (x — sv, v)dv 



L6/5 



T P" 



(4.19) 
(4.20) 



Here and in the sequel, the T -dependence of the constants C is continuous (on T G JR + ). 



Proof. The first estimate is (|4.17|l in Lemma 14.81 To derive the second one, we exploit eq. 
(14.181). the conservation of the L 2 -norm of the solution and (14.191): 



\Er(t)h < C s^(\\E^(t- S )\\ 2 + \\Er(t-s)h)\\w^(t-s)hds 
Jo 

< C||wo P || 2 sup sg(0ir] |s w Jwq p (x - sv,v)dv 
+ C\\w^\\ 2 f(t-8)-V 2 \\Er(8)\\2ds. 



1 

£2- 



The thesis follows by Gronwall's Lemma. 



□ 



<* _5 IK|| Ll(i 6/ 6) , t>0, (4.21) 



We shall now state a simple condition on wq that implies both conditions (A), (B). For 
wo € L\{L^]^) a Strichartz inequality for the free transport equation (cf. Th. 2 in |CP| ) 
reads: 

J wq(x — tv, v) dv 
and hence (A) and (B) hold. 

Remark 4.10 Let us again compare the a-priori bounds \^.iy\) , with their classical 

counterparts. Using j4-21\ ) we obtain the same t~ 2 -singularity of \\E wp (t)\\2 for the Wigner- 
Poisson system, as it was obtained in lOPf for the VP equation. In the latter case, similar 
L p -estimates hold for p in a non-trivial interval. One crucial reason for this difference is the 
conservation of L p -norm of the solution: while the WP equation only conserves the L 2 -norm, 
all LP-norms are constant in the VP case. A second reason is that we cannot exploit any 
pseudo-conformal law for the quantum case, since the Wigner functions are not non-negative 
(cf. \Pe\j for the classical case). 

As a by-product we obtain the following result for the self-consistent potential V, which 
follows directly from the splitting V" p = V wp + V{ p 



Vo P (x,t) 



Vr(x,t) 



/=i 1 1 

A E5s** W),0M). 



i=l 



(4.22) 
(4.23) 



Corollary 4.11 Under the assumptions of Proposition ^^ the self- consistent potential V n 
V Q wp + V™ p satisfies the following estimates for all t G (0, T] : 



IT7-UIP 



< cc T r 



< c 



l|2,SUp se(0jT] < S 



Wq p (x — sv, v)dv 



L6/5 



(4.24) 

,T\£- W . (4.25) 



20 



4.2 A-priori estimates for the electric field: the WPFP case 

According to Corollary 13.5) the mild solution of the WPFP problem satisfies for all t £ [0,T] 
(0 < T < t max ) 



w(x, V, t) 



+ 



G(t, x, v, xo,v )wo(x , v ) dx Q dv 



G(s, x, v, xq, vq)(@[V]w)(xo, vq, t — s) dxQ dvo ds 



with the Green's function G from 1)3.4 J) . According to jSCDMj we have 



f G(t,x,v,x ,v )dv = R(t)-^( X - X0 / Z ? i it)V0 

ir 3 V v-R(*j 



with 

M[x) 

m 

R(t) 



(2tt)" 3/2 exp 
1 - e-^ 



\x\ 
~2 



(3 

2at + a 



= 0(t), for t -» 0, 
4e _ # - e" 2/3i + - 3 



^ 3 



+ 4 7 



J 2 



(4.26) 
(4.27) 

0(t), fori^O. 

(4.28) 



By exploiting the redefinition 1)4.11)1 of the pseudo-differential operator, we obtain the fol- 
lowing expression for the density n[w] 



n[w] (x, t) 



w(x, v, t) dv 



1 



R{tfl 2 
rt i 



X — Xq — d(t)vQ 



w (xo,v Q )dx dvo 



+ 



R(s) 3 / 2 



X — Xq — $(s)vq 



yfRjs) 



div VQ (T[V X0 V]w) (x ,v ,t- s) dx dv ds 



1 



( X - X \ 



R(t)WJJ \VW)J 



+ 



w (x - tf(t)v , v ) dx dv 

X — Xq — 1?(s)Uq 



o R(s) 2 
n (x,t) + ni(x,t), 



(T[V X0 V]w) (x ,v ,t- s) dx Q dvo ds 



where 
n (x,t) :-- 

ni(x,t) :-- 



i?(s) 



R{sf/ 2 \^/Rjs~ 
Correspondingly, we can split the field (with A = ^-): 

E (x,t) := *x n (x,t) 



* x riQ (x, t) , with ng(x,t):= / wq{x — ~d(t)v, v) dv , 
* x div x / ^rfV^yjiuJ (x — i?(s)u, v,t — s) dv ds. 



x 



R{tf/ 2 \^W) 



* x Eg(x,t), (4.29) 



21 



with 



E$(x,t) := A-^r * x no(x,t), 



Ei(x,t) := \- T r* x ni(x,t). 



(4.30) 



Remark 4.12 Note that the splitting of the density (and of the electric field) is the same 
as in Wo 11 \Bo2\\Ga2^ : in the WPFP case the two components of the decomposition (no,m, 
as well as Eq,E±) are smoothed versions (in fact, convolutions with a Gaussian) of those 
appearing in the WP case (namely n^ p , n± p , E^ p , E™ p ) . Actually, the density tiq(x ,t) (which is 
convoluted with the Gaussian to give no) already differs from Hq p {x, t) in the WP case because 
the shift contains the function which is due to friction (and analogously for EQ(x,t) and 
E^(x,t)). 

From Lemma 14.61 we directly get 
(r[i?]n) (x — i}(s)v, v, t — s) dv 



Ml 



La(K|) 



< C${s)- z l 2 \\E(t-s)\\ 2 \\u(t-s)\\2, Vt>s>0. 

(4.31) 

To derive an L 2 -estimate on the field we shall proceed as in the WP case (Lemma 14.81 
Proposition 14.9)1 . 

Lemma 4.13 Let w be the mild solution of the WPFP equation \3.7j) and let wq £ X satisfy 
JX]) for some uj £ [0, 1). For any fixed T > the electric field then satisfies V X V £ V Tu) _i 
and the following estimates hold: 

1. for2<p< 6, 



30-2) 
2p 



\E (t)\\p < C^lkolllll^Wliy^Or-M), Vt 6(0,3*]; 



(4.32) 



2. 



\E, 



2 < C [T, \\wq\\ 2 , sup 

V 86(0,71 



sr\\4(s)\\ L6/5 



})**" 



Vt£(0,T]. (4.33) 



Proof. The estimate for ||i?o(i)||p, p £ [2,6] is obtained by applying first the generalized 
Young inequality and then the Young inequality to the expression (|4.29j) 



\\Eo(t)\\ P < C 
< C 



1 



R(t)y 
i 



X 



* x no(x,t) 



\4( x ,t)\\ c 



CM 



'/• 



with q 



3p 



£ [6/5,2]. 



p + 3 

Next we interpolate between I? and L 6//5 , use Ij2.3j) and the dissipativity of the operator 

3 
2 



-v V x - | in X (cf. Lemma l2~Sl): 



11 r 



with 



f — ~. Hence 

z q 



< C\\w o (x-0(t)v,v)\\ e x \\T$ 

< Cel^\\w \\Unt(t)\\lj 5 6 , 



\E (t)\\ p < c(r)|KI&ll4(*)fe/i 



1-0 

6/5 
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We rewrite the function Ei(x,t) as 



k=l 



■AT 



F k (x,t,s)ds, (4.34) 



o R(s) 3 / 2 \VW)J 

with Fk(x,t,s) := / ' (Tk[Eo + E±]w) (x — $(s)v,v,t — s) dv, 

For estimating it we exploit classical properties of the convolution with the kernel ^ and 
apply the Young inequality: 



\Ei(t)\\ 2 < C U{a) 



< C / 



< C(T)\\w \\ 2 



* \\E (t - s)\\ 2 + ||^i (t-s)|| 2 



* x F(x,t,s) ds 

2 

F(x,t,s)\\2ds 
ds, 



where the last inequality follows from (|4.31l) and the L 2 -a-priori estimate on the solution w 
(cf. Lemma l4~2|) . By applying the estimate (|4.32[) to ||£?o(i)||2 ; we get 



|£i(*)||2 < C(T)|K||2 sup 

\te(o,T] 

< \\Ei(t ~ s)h 



tlr 



L 6/b \ / 0(s)-30(t-«)- w ds 



+ 



v 7 ^) 



ds 



(4.35) 



where the function i9(s) = 0{s) as s — > 0. Thus the integrals are finite. 
To establish a solution of (|4.34j) we introduce the fixed point map 



(ME)j(x,t) := A£ 



1 ^ M' x 



*x / (Tk[E() + E]w) (x — $(s)v, v,t — s) dv ds. 



By using < ^p-t < $(t), Vi G (0, T] and (|4,35|) . a simple fixed point argument as in the 
proof of Lemmma 14.81 with the contractivity estimate: 

7rfr(|-o;) 



\M n E(t)-M n E(t)\\ 2 < [C ] 



n + l 



\woh] t 2 



sup (s^\\E(s) - E(s)\\ 2 ) 



shows that the linear equation (|4.34|) has a unique solution E\ G V T i. Hence V X V = 
Eq + Ei G V T w _ i and Gronwall's Lemma then yields estimate (|4.33|) . rj 



Remark 4.14 For the derivation of the a-priori bound on\\E\\2, we did not use any moments 
ofw (neither inx norv), nor pseudo-conformal laws (cf. \Bol\ \Bo^ \Pe\ HJaty for the classical 
analogue, i.e. VPFP). In fact, the latter are not useful in the quantum case, since the Wigner 
function typically also takes negative values. Moreover, the convolution with the Gaussian 
did not play a role there; the estimate \4-^'^ relies just on the dispersive effect of the free- 
streaming operator. The parabolic regularization will be exploited in the "post-processing" 
Proposition \4-15\ 
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The above lemma was the first crucial step towards proving global existence of the WPFP 
solution. Next we shall extend this estimates on the field £ to a range of L p -norms: 

Proposition 4.15 Let w be the mild solution of the WPFP equation {3. 7| ) and let wq € X 

satisfy iTXj) for some uj £ [0, 1). Then, we have for any fixed T > and for all p G [2, 6): 



< c r.iKlb, sup {Hsr\\nt( s )\\ L6/5 })t- P 



3_ I 
4 " 



Vt€(0,T]. (4.36) 



PROOF. We shall estimate Ei(t) (cf. (E34m by using classical properties of the convolution 
by the kernel ^ and the following 



CR(s)^5, VI < q < oo. 



(4.37) 



Namely, 

l|£iM(t)ILp(H») < c/V^) 
< c 



i?( S ) 3 / 2 iv^J 



* x F(x,t,s) 



LP (Ml) 



ds 



\Eo(t - s)\\ 2 + ||^iH(* - s)||2) IH* - s)\\ 2 ds, 



VW) 

where we used the Young inequality with 1 + 1/p = 1/q + 1/2 (thus, p > 2) and, for the 
L 2 -norm, the Lemma l4.(il Then, by applying Lemma 14.131 we get 

II -El H(*)IIlp(]r3) < c(t, su^{i?(*) w ||n^(s)|| L 6/5} J|wo||i^(]R»)V 



VW) 



t - s)-" + (t - ds 



Since 0(t) = 0(t),R(t) = 0{t) for t -> (cf. 023, @35» , the last integral is finite for 
all t > and for 3/(2g) - 2 > -1 44> 3/(2p) - 5/4 > -1 44> p < 6. In fact the integral is 
0(ti-^). D 

Remark 4.16 Proposition \4-l$ provides a non-trivial interval of LP- estimates for the elec- 
tric field in the WPFP case. This is due to the regularizing effect of the FP term. We 
remark that the corresponding Gaussian is "better behaved" than the classical one, since the 
quantum FP operator is uniformly elliptic in both x and v variables. On the other hand, 
exactly as in the WP case, the range of LP- estimates for the WPFP equation is smaller in 
comparison to the counterpart VPFP and that depends again on the non- negativity of the 
classical distribution function. 

As a further result, we obtain an a-posteriori information on the self-consistent potential 
V, which follows directly from the a-priori estimates on the field. Accordingly, we split the 
potential as V = Vq + V±, with 

3 

V (x,t) := XY / ^h*AE ) i (x,t), (4.38) 
i=i 1 1 
3 

V r {x,t) := A^^^^t). (4.39) 
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Corollary 4.17 Under the assumptions of Proposition \4-15\ the self- consistent potential 
V(t) = V (t) + Vi(t) belongs to L P (JR 3 X ) with 6 < p < oo, and satisfies for all t G (0, T] : 

\\V (t)\\p < C(^)ll^o||Scl«(*)llia; 5 = with6 = \--, (4.40) 

\\Vi(t)\\ P < C(r,|K|| 2> sup {^ni4( S )|| i6/5 }^i + ^ w . (4.41) 

4.3 A-priori estimates for the weighted L 2 -norms 

A first consequence of the a-priori estimates for the electric field is the following 

Lemma 4.18 For all wq E X such that IX)) holds for some u> £ [0, 1), the mild solution of 
the WPFP equation \3. 7| ) satisfies 



-fill 



< C[T,\\w \\x, sup U(sr\\4(s)\\ L(i/5 \ ) , Vie [0,21. (4.42) 



Proof. In order to justify the derivation of this a-priori estimate we need again the approx- 
imating classical solutions y n introduced in the proof of Lemma 14.21 Mutiplying both sides 
of ()4.3|) by v 2 y n (t) and integrating yields 

1 d " —12 f f 2 / + wr 



2^lMnWI| 2 = JJv^n^Ayn^dxdv + JJviyn^fn^dxdv. 
By analogous calculations as in the proof of Lemma 12.81 (cf. also (j5.3|) ) we get, 

R 

\v\ 2 y n (t)Ay n (t)dxdv < 3a\\y n (t)\\l + -\\vy n (t)\\l, 
and hence 

^Hn(t)||i < M\Vn{t)\\l + ^\\vy n {t)\\l + JJ\v\ 2 y n {t)f n {t)dxdv, Vte[0,T]. 

By integrating in t, letting n — > oo, and using (|4.1|) . we have 

2(T /"* 

< ||^o||| + -g(e 3/3 *-l)IKI|| + /9/JMs)|||d S 



o 



ft 

|2 



+ 2y JJ \v\ 2 w(s)f(s)dxdvds, Vie[0,T]. 

Using again the skew-symmetry of the pseudo-differential operator and the Holder inequality 
yields 

J JJ Viw(s)vif(s) dx dv ds = -J JJ Viw(s)U[diV[w(s)]]w(s) dx dv ds 

< I [\\ViW(8)\\2\\n[aiV[w(8)]]w(8)\\2ds, 

* Jo 

with the operator £1 defined in (j2.10j) . Estimating as in l|2.12j) and using the Sobolev inequal- 
ity we obtain for t £ [0, T]: 

\\mviw(t)]] W (t)\\ L2(JR 3 XjR3v) < cpiViwmwmLHKixjRS) 

< C\\aiV[w(t)]]h\\vv(jt)h, (4-43) 
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where w(x,r),t) := J : v ^ v (w(x,v,t)). Finally, using Prop. 14.151 (estimate 1)4.36)1 with p = 3) 
yields 

\\vw(t)\\ 2 2 < C(T) (llvwolll + H^olll) +C H^oll^, sup {0(*) w ||ng(«) \\ L e/ B fj 

'(a - * + s -«>+i +p\ \\vw(s) ||| da, te[0,T], (4.44) 
and the Gronwall Lemma gives the result. rj 

With this result we can proceed to derive the a-priori estimate for |||w| 2 u)(t)||2. 

Lemma 4.19 For all wq G X such that IX)) holds for some to G [0, 1), i/ie miZcZ solution of 
the WPFP equation j3.7\) satisfies 

\\\v\ 2 w(t)\\ 2 2 < c(t,\\w \\ x , sup U(sr\\ni(s)\\ L6/5 \) , VtG[0,T]. (4.45) 

Proof. In order to control the term |||v| 2 u;(i)||2, we shall use the same strategy as in the 
Lemmata 14. 21 and 14. 181 Multiplying both sides of 1)4.3 J) by vfy n (t) and integrating we get by 
using ()5.3|) and repeating the same limit procedure as in the previous lemma: 

\^hjjviv{t?dxdv < Vo\\wm\l + (zo- l -^Y.jjvtw{tfdxdv 

3 

+ J2[[vtw(t)f(t)dxdv, VtG[0,T]. 

1=1 

By integrating in t, using Ci|t>| 4 < J2 v t ^ C2 1 ^ | and 1)4.1)) . we have 



\v\Mt)\\l 



< c(\\\v\ 2 w \\ 2 2 + j(e 3f3t - 1)|KH1 + (6a - P)j\\\v\ 2 w(s)\\l ds 

+2Y,f t ffvtw(s)f(s)dxdvds], Vte[0,T]. (4.46) 



3 „j 

Using again the skew-symmetry of the pseudo-differential operator 0, the equation 1)2.9)1 and 
the Holder inequality, we have 

v 2 w{s)v 2 f{s)dxdvds < \ f^vfwis^WeidfViwis^w^hds 

4Jo 

+ / \\vMs)h MdiVHs^w^hds. (4.47) 
Jo 

Since u>(x, . ,t) G i? 2 (K^), the Gagliardo-Nirenberg inequality yields for t G [0, T] 



\w(x, . ,t)\\ L oo {m 3) < C\\w(x, . ,*)Hi6 (JR 3 ) || \v\ 2 w(x, . ,t)||^2 2 (]R 3 



< Clltw^.^Jll^jllHawCx,.,*)!^ 2 ^. (4.48) 



Using 

\\&V[w(t))\\ 2 = \\n[w(t))\\ 2 = C||^(.,r / = 0,t)|| L2(]R 3 ) < C^J\\w(x, . , i)!!^^ efc 



1/2 
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(|4.48|) . the Holder inequality, and (|4.42() we can estimate: 

v 1/2 



\@[dfV[w(t)]]w(t)\\ 2 < C||AyKi)]|| 2 ^||«i( a ;,.,t)||^ 00(]R s ) dx^ 
< C j\\w{x,.,t)\\lodx 

\VW(X, • ,^1^2(^3) \\\V\ 2 W(X, . ,01^2(^3) dx 



C 



< C[T,\\wo\\x, sup {^sr\\4(s)\\ L6/5 })\\\v\ 2 w(t)h- (4.49) 
For the second term of the r.h.s. of (|4.47|) we proceed as in (j4.43|) and use the estimate 

WnidiVMmviwWh < c\\a i v[w(t)]\\ 3 \\mmv'(^ i ifi(B')) 



< C\T,\\w \\x, sup 

V 86(0,71 



{^(«ri|no(*)llL«/«}J (t^+t~^) \\\v\ 2 w(t)h, Vte[0,T\. 



(4.50) 

Analogously to Q4.44|) , combining the estimates ()4.47j) , (|4.49|) and Q4.5U|) the Gronwall Lemma 
gives the assertion. rj 

Proof of Theorem 14.11 

The Lemmata 14.21 and 14. 191 show that 



\w(t)\\ x < C[T,\\w \\x, sup U(sr\\nt(s)\\ L6/5 \ , We[0,T], V0 < T < t 
V se(o,T] J J 



max ■ 



with C being continuous in T £ [0,t max ]. Then, Corollary 13.51 shows that the mild solution 
w exists on [0, oo). rj 



4.4 Regularity 

The following result concerns the smoothness of the solution of WPFP, the macroscopic 
density and the force field, for positive times. 



Corollary 4.20 Under the assumptions of Theorem \4-l\ the mild solution of the WPFP 
equation \3.7\) satisfies 

weC((0,oo);C£>(]R 6 )), 

n(t),E(t),V(t) eC((0,oc);C^(IR 3 )). 

Proof. Obviously, w(t) £ C(IR 6 ) Vi > 0, because of the Green's function representation in 
(|3.22j) . (|3.3|) . If we differentiate equation (|3.7|) with respect to xi and, resp., Vi, we obtain 
the following linear, inhomogeneous problems for any fixed t\ > 0. 

z t (t) = Az(t) + e[V[z(t)]]w(t) + &[V[w(t)]]z(t), Vt>ti, z(h) = d Xi w(h) G X, 
y t {t) = Ay(t)+f3y(t)-d x Mt) + ®[V[w(t)}]y{t), Vt > *i, y(h) = d v .w(h) € X. 



27 



By arguments analogous to Lemma 13,31 there exists a unique mild solution 

z = 8 Xi w G C{[t 1 ,oo);H 1 (M 6 ; (1 + \v\ 2 ) 2 dxdv)). (4.51) 

By an induction procedure, the derivatives V^V^w, for a, (3 G IN 3 , |a| + |/?| = m > 1 
are also mild solutions of similar problems with additional well-defined inhomogeneities and 
with initial times < t\ < t 2 < ... < t m . This yields V£V^w G C([i m , oo); fl^IR 6 ; (1 + 
|f| 2 ) 2 eZt>)), and thus V"V^w G C((0, oo);X). Hence, the statement about smoothness of 
the density and the electric field is straightforward from Propositions l2.ll and l2.3l and Sobolev 
embeddings. n 



5 Appendix 
Proof of Lemma 12.81 

For u G D(A) we have 



< Au, u >x = < Au,u >i2( K 6) +^ vfuAu 



i=l 



(5.1) 



where /// denotes the integral Jjps J^f^x, v) dv dx , and the norm ||. ||^ is defined by (|2.7j) . 
Using integrations by parts we shall calculate the three terms on the right hand side sepa- 
rately. 



< Au,u > L 2 (]R 6' ) = ^ ( - ViU Xi u + (3 (viu) Vi u 



3 



< E 

i=i 



3(3 II u 2 + (3 ViU Vi u -ail u 2 



+7 e + 



»7 



a //u 



— >0||xft||i + Q - a ) ||V„ti||! + (£7 - a) IIV^IH- 



With e = r we obtain 



< Au,u > L2(]R 6 ) < 
Next we estimate the second term of (|5.1[) : 

Y.ff v i uAu = E (-JJvtv j u X] u + f3jJvt(v j u) Vj u 

i=l i,j=l 

+cr vfu v „ v „u + 27 vfu XjV ,u + a vfu XiXj u 



< E 

4 



(3 v { u +(3 ViVjUv.u-a vtu v 



-a Uvfu Vi u + 7 e //t;*«L + 7 



a I Ivfu 2 



(5.2) 
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3 , ! 



1=1 



< 9ff||«||l+ f-^ + 3a)W/^ 2 , (5.3) 



i=i 



by choosing e = - and by an interpolation. 
Collecting the two estimates yields 



<Au,u >jc < Q/3 + 9a) \\u\\l + 3a YjjvW < + 9^ \\ufg. 

Thus, the operator ^4 — kI is dissipative. □ 

Proof of Lemma 12.101 

To prove the assertion we shall construct for each / G D(P) C L 2 (IR 6 ) a sequence {/ n } C 
D(P) such that / n — > / in the graph norm 

II/IIp = ll/lk* + IIM 2 /IIl 2 + IIM 2 /IIl 2 + IIP/IIl" + IIM^/II^ + IIM 2 ^/IIl- 

To shorten the proof we shall consider here only the case 

P = e + uv ■ V x + \tx ■ V„ + (3v ■ V„ + aAj + aA t , + jdiv v V x 

(cf. the definition of the operator A in (|2.13|l ). but exactly the same strategy extends to the 
case, where P is a general quadratic polynomial. 
First we define the mollifying delta sequence 

4> n (x,v) : = n 6 4>{nx , nv) , n G IN, G IR 3 , 

where 

0GC o °°(IR 6 ), 0(a;,t;)>O, 

(x,v)dxdv = 1, and supp0 C {\x\ 2 + |t>| 2 < 1}. 



By definition we have the following properties: 
(I) ^^inP'OR 6 ), 

(II) ld Xi 4> n , ld Vi cj> n ^ in P'(IR 6 ), i = 1,2, 3, 

(III) (x, u) a fl2. v) [(x,u) 7 n (a;,v)] ^ in £>'(IR 6 ), with a,/?, 7 G 1N§ multi-indexes 

and |7| > 0, since (x,i;) 7 n -> in £>'(IR 6 ). 
The cutoff sequence is 

i/> n (x, v) := V ( ~~ ) ' n G IN, x,v £ IR 3 , 

where tp satisfies 

y> g c °°(ir), < ^(*) < 1, supp^ c [-1, 1], = 1, 

and 

\^(z)\<Cj, VzGM, j = 1,2. 
The sequence ^ n has the following properties: 
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(IV) ip n — > 1 pointwise, 

(V) (z,v)"tf x ^(z,v) = ^^gyV (^) , with a,/3 € INg, 
l a l = l/^l = 1) are supported in the annulus 

supp(V(3?)) = {{x,v)\n/2<\(x,v)\<n} =: V n , 
and they are in L°°(]R 6 ), uniformly in n G IN. 
(VI) nd^nO^) = jg^' (M) , with a E ]Ng, \a\ = 1, 
are uniformly bounded in L°°(IR 6 ). 

(VII) ^^VnC^.w) = J^V^" (^T^) + (n y |(lc,ti)| ~ Jk^)i 3 ) ^' (^^) ' 
with |a| = 2 have support on V n and converge uniformly to in L°°(1R 6 

We now define the approximating sequence 

fn{x,v) := (f * (f> n )(x,v) ■ 1p n (x,v), 11 £ IN, 

where '*' denotes the convolution in x and v. 
By construction we have f n £ Co°(IR 6 ) = D(P). 
Since we can split our operator as 

3 r0 

- + uvid x . + nxid Vi + (3vid v . + ad 2 . + ad\ + jd v .d Xi 



E 

8=1 

3 



i=l 

we shall in the sequel only consider 

P = p{y,z,d y ,d z ), y,z e IR 

acting in one spatial direction y = Xj and one velocity direction z = Vj. 
We have to prove that f n (x, v) — ► /(x, u) in the graph norm 



L2 + HM 2 /lb + IIM 2 /lb + WPfh* + \\\A 2 Pf\\L* + \\\v?Pf\\v- 

According to the 6 terms of the graph norm we split the proof into 6 steps: 
Step 1: By applying (PI) and (P4), we have 

/„-/ in L 2 (K 6 ). 



Step 2: For the second term of the graph norm we write 

xjfn = {xlf * <j) n )lpn + 2(xj/ * Xi4> n )lp n + (/ * xf<j) n )ljj n . 

The first summand converges to xff in L 2 (IR 6 ) and both the second and the third terms 
converge to by (III), since also Xif belongs to L 2 (IR 6 ) by interpolation. 

Step 3: For the third term of the graph norm the same argument as in previous step can 
be used. Hence we have 

f n -> / in Z. 
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Step 4: To prove that Pf n -> Pf in L 2 (IR 6 ) we write: 


Pfn = -(/ * 4>n)lpn + "(zfy * <Pn)lpn + Kvfz* fintyn + P(zfz* <Pn)lpn 
+a(f yy * <t>n)i>n + °{fzz * 4>n)?Pn + l{fyz * <j>n)i>n + r„(y, z) 

= (Pf*^ n )i> n +r l n {y,z). 

As we shall show, all thirteen terms of the remainder 

r n = v{f * d y (z(j) n ))tp n + u(f * (j) n )zd y ilj n + fi(f * yd z (p n )ijj n 

+fi(f * 4> n )yd z ^n + 0(f * d z {z^ n ))lp n + 0(f * 4> n ))zd z ll) n 

+2a(f * (-dy4> n ))(ndyi(; n ) + a(f * 4> n ))(d 2 ip n ) + 2a (f * -d z (f) n )nd z ip n 
n y n 

+<r(f * 4>n)d z ^ n + t(/ * (-d z 4>n))(nd y i) n ) + j(f * (-d y (j) n )){nd z i) n ) 

n n 

+l{f * K)dyd z ^ n 

converge to in L 2 (IR 6 ). 

The first, the third and the fifth terms converge to in L 2 (1R 6 ) by (III). 
In the second, fourth and the sixth terms, exploiting (V) we have 

||(/ * cP n )(zdy^ n )\\ L 2 {R6) < C\\f * <j> n - f\\ L 2 {Vn) + ||/|U 2(yn) - 0, (5.4) 
because ||/|| L 2(R6) = ||/||l2(b 1/2 (o)) + SZfcLo ll/IU 2 (v^ fc )- 

For what the seventh, ninth, eleventh and twelfth terms are concerned, we can exploit (VI) 
and then (II). 

The remaining terms can be handled thanks to (VII). 

Step 5: To prove that \x\ 2 Pf n —* \x\ 2 Pf in L 2 (IR 6 ) we write: 

Q 

xfPfn = - {xff * <p n )lp n + v(x\zf y * <f> n )lp n + H(x 2 yf z * 4>n)^n + P(x 2 zf z * 4> n )^n 

+a(x 2 f yy * <t> n )ip n + a(x 2 f zz * (j} n )ip n + l{x 2 f yz * <j> n )i])n + rl(y, z) 
= ( x fPf *4>n)ipn + r 2 (y,z). 

The remainder r 2 can be split in the following way (y = Xj, z = Vj): 
2 6 

r 2 n v = 2u(zxif * d y (xi(J) n ))ip n - 2u8ij(zf * Xi<f> n )ip n + u(zf * dy{x 2 4> n ))il) n 
+v{x 2 f * zd y (fi n )ip n + 2v(xif * Xizdy(f) n )i/) n + u(f * x 2 zd y <p n )ip n 
+v{x 2 f * 4> n )zd y ip n + 2v(xif * Xi<t> n )zd y ip n + v(f * x 2 <j) n ) zd y i/j n 

r \,a = 2/z(x i y/*^(x i ^ n ,))V'„ + /i(y/*9 z (x 2 (/> n ))V'n + M(^ 2 /*y^0 n )V'n 
+2fi(xif * Xiyd z cp n )ip n + n(f * x 2 yd z 4> n )^ n + fi(x 2 f * <pn)yd z ip n 
+2fi(xif * Xi(f) n )yd z ip n + n(f * x 2 4> n )yd z ip n 

r 2 j/3 = 2f3(xizf * Xid z (f> n )ip n - 2(5(xif * Xi(p n )ip n + (3(zf * x 2 d z 4> n )ip n - f3(f * x 2 cp n )ip n 
+/3(x 2 f * d z (z4> n ))i) n + 2/3(xif * Xid z (z$ n ))i) n + /3(f * x 2 d z (z<f) n ))i) n 
+(3(x 2 f * 4> n )zd z i) n + 2(3(xif * Xi4> n )zd z ^ n + 0(f * x 2 4> n ) zd z i/j n 
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r 2 na = 2a(xif * dy(xi(j) n ))ip n - 4a<%(/ * d y (xi4) n ))^ n + a(f * dy(x^ n )^ n 

+2a(x 2 f * -d y (j) n )ndyij} n + 4a(xif * -^d y (j) n )ndyip n + 2a(f * -^d y 4> n )nd y i) n 

+a(x 2 f * 4>n)d y tp n + 2a(xif * Xi(j) n )d^ n + a(f * x 2 (j) n )d 2 ip n 

r 2 nry = 2a{xif * Xid 2 (f) n )4} n + a(f * x 2 d 2 z (f) n )ip n + 2a(x 2 f * -d z 4> n )nd z i) n 

n 

x x 2 
+Aa(xif * —d z 4> n )nd z il) n + 2a(f * —d z (/) n )nd z ijj n + a(x 2 f * § n )d 2 z ^ n 
n n 

+2a(xif * Xi(p n )d 2 ip n + a(f * x 2 (j) n )d 2 ijj n 
r\ a = 2^{xif * d y (xid z (j) n ))ip n - 27<%(/ * Xid z (j) n )ij) n + j(f * d y d z (x 2 (j) n ))ip n 

+l( x if * -d z (j> n )ndyil) n + 2>y(xif * —d z (p n )nd y ip n + -y(f * —d z (j) n )ndyij; n 
n n n 

+l( x if * -dy(t)n)nd z if; n + 2j(xif * —d y $ n )nd z i\) n + -y(f * —d y (f> n )nd z ^ n 
n n n 

+l( x lf * 4>n)dyd z ip n + 2j(xif * Xi4> n )d y d z ip n + 7(/ * x 2 4> n )d y d z ip n . 

By the properties (I)-(VII) and estimate like (|5,4|) . it can be easily seen that each term con- 
verges to in L 2 (R. 6 ). 

Step 6: In analogy to \x\ 2 Pf n , the sequence \v\ 2 Pf can be split as 

v 2 Pf n = (v 2 Pf * K)^n + r 3 n (y, z). 

Due to the symmetry of the operator P in x and v , the terms of the remainder r 3 can be 
obtained from r\ by interchanging y and z (and changing the coefficients), except for the 
following term 

vf[Pzd z ((f *(j>n)lpn)] = P(v 2 zf z *<j) n )iJ n + r^ f3 , 

where 

r z n f3 = 2/3(vizf * d z (vi<t> n ))il) n - 2/3(1 + Sij){vif * ^„)V>n + @{zf * d z {v 2 <t) n ))ipn 
-f3(f * v 2 (f) n )^ n + (3(v 2 f * d z (zcf) n ))4) n + 2(3(vif * Vid z {z4> n ))i) n 
+/?(/ * v 2 d z (z(/) n ))^ n + (5{v 2 f * 4> n )zd z ij} n + 2f3(vif * Vi(f) n )zd z ip n 
+/?(/ * v 2 4> n )zd z ^ n 

converges to in L 2 (IR 6 ), since (I)-(VII) and (|5.4|) can be used. □ 

Proof of Proposition 13.11 

First, we shall prove the following estimates on the derivatives of the Green's function (|8.4jl : 

C(f - - 2& ^2.) 

\V v G(t,x,v,x ,v )\ < b 2 '^ 2 ' 2 , Vt<i , (5.5) 

C(f % v XQ_ vo\ 

\V x G(t,x,v,x ,v )\ < b' l ' 2 '^ 2 ' 2 ~ , Vt<h. (5.6) 

with b = 6(q, 7, a), to = to(a, (3, a, 7), b' = b'(a, 7, a) and t\ = t±(a, (3, a, 7). The ^-derivative 
of G is given by 



V v G(t,x,v,x , vq) = G(t,x,v,x ,v ) 



^(t)e^ - 2v{t)^f) (x - - x ) 



/(*) 
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(5.7) 



2A(t)e# - ^^f 1 ) (e&v - v ) 

W) 

For all real a, b, c > such that c/\/a < by2e, one easily verifies that 

c\x\ < be a ^ 2 , VxeK 3 . 
Since a, a > 0, we have for t > small enough 

KO-^W > °> A(t)--Mt) > o. 

In order to apply the estimate (|5.8() to the two terms inside the squared bracket in 1)5.7(1 we 
shall use for t small: 



(5.8) 



Cl 


Vt 

m 


H(t)ef* - 






2 7 








/3"(*)- 
V4 / 


-in® 
it) 




-7 2 )(a 


+ 7) 


= ll y/3( acr ~ 7 2 )( cr + 7) 


Similarly, 








C2 


Vt 

m 


2A(t)e /3t - 


-Kty-^ 




2a 




\f02 






\^t) 


y / 3(aa 


- 7 2 )(a 


+ 7) 



< &i\/2e, 



4 /(*) 



< 6 2 V2^, 



with 62 = a/^3(aa — 7 2 )(a + 7). Then, there exists some to > such that, for all t < to, 
the two inequalities can be combined with b = max{6i, 62} to give 



L(t)e^ - 2v(t)'- 



( x _ e^-l 



v-x ) + (2A(t)e^ - li{t)^) (ePv - v ) 



f(t) 



Vt < 



< 6exp 



H(t)e^ - 2u(t) 







x — V — XQ 



+ 



2A(i)e /3i - n{t) 



+ (\{t) - \^t)) 



if® 



< b 



u{t) 



exp 



V -Xo 



+ A(t) 



e^u — Uq 



+ 



s/(*) 



le^v - vo\ 



v — xq I • I e^*f — vq 



Hence, 



\V v G(t,x,v,x ,v )\ < b 



G(t,X, V,Xq,Vq) 

Vt 



ex P \ 7" 



v{t) 


X — V - X 


2 + A(t) 


e^*t> — 


2 + /x(t) (x - ^^-v - JC ) 


■ 




f(t) 



and the decay (|5.5|) follows by comparison with 1)3.4)1 . 
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Next we consider the x-derivative of the Green's function, 
V x G(t,x,v,x ,v ) = G(t,x,v,x ,v 



2v(t){x - ^g^v - x ) + fi(t)(e^v - vo) 



f(t) 



Analogously, the decay 1)5. 6 j) follows by exploiting that for t small enough 

2a 



4 



^/3(aa -7 2 )( CJ + 7) 



2;. 



3 x(t)-y(t) ^3(aa -7 2 )(« + 7) 



< 6' 2 V2e, 



m 

with appropriate b'i(a, 7, a), b' 2 (a, 7, a) 
Since 



i A / \ 

e ruo(a;,v) 



G(t, x, v, x ,vo)wo(x , v ) dx dv , 



we have 



V v e w (x,v)\ < J J I V v G(t,x,v,x , v )\ \w (x , v )\ dx dv 

7 -1/2 ff^ ( x v x v o\ 1 / \ 1 

< bt 1 G \t, j \wo{x ,v )\dx dv 



64bt ^ 2 J J G(t, x, v, x~o, vo)\wo{2x~o, 2vq)\ dxo dvQ 
64br 1/2 e tA w (x,v), Vi < t . 



(5.9) 



Here we used the decay (|5.5j) . and we put x = §, u = | and u)o(x,v) = |u>o(2x, 2v)|. 
The assertion (|3.5f) follows directly by applying the estimate ()3.2|) to ()5.9|) and choosing 
T = min{t ,ti}- 

The estimate (|3,6(1 can be obtained analogously. □ 
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